Alternative Methods to Finding Patterns in 
HiRes Stereo Data. 



R.U. Abbasi/'* T. Abu-Zayyad/ J.F. Amman, ^ 
G.C. Archbold, ^ K. Belov, S.A. Blake, ^ J.W. Belz, ^ 

O : S. BenZvi, ^ D.R. Bergman, ^ J.H. Boyer, ^ G.W. Burt, ^ 

^ : Z. Cao, ^ B.M. Connolly, ^ W. Deng, ^ Y. Fedorova, ^ J. Findlay, ^ 

X) : C.B. Finley, ^ R.C. Gray ^ W.F. Hanlon, ^ CM. Hoffman, ^ 

[T^ ■ M.H. Holzscheiter, ^ G.A. Hughes, P. Hiintemeyer, ^ 

m ; B.F Jones ^ C.C.H. Jul, ^ K. Kim, ^ M. A. Kirn, ^ B.C. Knapp, ^ 

^ ; E.G. Loh ^ M.M. Maestas, ^ N. Manago, ^ E.J. Mannel, ^ 

^ : L.J. Marek, ^ K. Martins, ^ J.A.J. Matthews, ^ J.N. Matthews, ^ 

^ : S.A. Moore ^ A. O'Neill, ^ G.A. Painter, L. Perera, ^ K. Reil, ^ 

m : R. Riehle, ^ M. Roberts, § D. Rodriguez, ^ M. Sasaki, ^ 

O ! S. Schnetzer,"^ L.M. Scott, ^ M. Seman, G. Sinnis, 

O ■ J-D- Smith, P. Sokolsky, G. Song, ^ R.W. Springer, 

' B.T. Stokes, ^ J.R. Thomas, ^ S.B. Thomas, ^ G.B. Thomson, ^ 

O ■ 

' '^University of Utah, Department of Physics and High Energy Astrophysics 

Ci '■ Institute, Salt Lake City, UT 84112, USA 

^ ■ ^Los Alamos National Laboratory, Los Alamos, NM 87545, USA 

^ ' ^ University of Adelaide, Department of Physics, Adelaide, South Australia, 

Australia 

'^Montana State University, Department of Physics, Bozeman, Montana, USA 

^The State University of New Jersey, Department of Physics and Astronomy, 

Piscataway, New Jersey, USA 

^ Columbia University, Department of Physics and Nevis Laboratory, New York, 

New York, USA 

^University of New Mexico, Department of Physics and Astronomy, Albuquerque, 

New Mexico, USA 

^University of Tokyo, Institute for Cosmic Ray Research, Kashiwa, Japan 



D. Tupa,'' S. Westerhoff/ L.R. Weincke," A. Zech," X. Zhang,' 



Preprint submitted to Elsevier 



5 February 2008 



Abstract 



In this paper Ultra High Energy Cosmic Rays UHECRs data observed by the HiRes 
fluorescence detector in stereo mode is analyzed to search for events in the sky with 
an arrival direction lying on a great circle. Such structure is known as the arc 
structure. The arc structure is expected when the charged cosmic rays pass through 
the galactic magnetic field. The arcs searched for could represent a broad or a small 
scale anisotropy depending on the proposed source model for the UHECRs. The 
Arcs in this paper are looked for using Hough transform were Hough transform is 
a technique used to looking for patterns in images. No statistically significant arcs 
were found in this study. 

Key words: Anisotropy, Hough Transform, Ultra high-energy cosmic rays, HiRes, 
AGASA 



1 Introduction 

Cosmic rays passing through the galactic magnetic field are deflected. The 
degree of deflection depends upon their energy and the strength and structure 
of the field. This deflection has the potential of producing either large scale 
or small scale anisotropy in the cosmic ray arrival direction. The scale of this 
anisotropy is determined by the UHECR source model. In particular, with 
moderate deflection through a coherent field, the pointing directions of the 
deflected cosmic rays from a point source may form an arc in the sky, as 
illustrated in Figure 1 [1]. Alternately, an extended candidate source, such as 
the Galaxy, may lead to a lengthened circle of events. These are the type of 
structures we would like to identify. 

In this study we use the Hough transform to search for arc-like structure. 
Hough transform is a class of duality transformations commonly used to iden- 
tify 1-dimensional structures in 2-dimensional graphical images [2] . The goal 
of this study is to look for arcs in the HiRes stereo data set. Here we use the 
term "arc" to describe a segment of great circles across the sphere of the sky. 

In this study, the Hough transform is applied as a self- dual mapping on a 
sphere, where the pointing direction of each event is mapped to a great circle 
that is perpendicular to its arrival direction. For events that lie on an arc, their 
associated great circles would intersect at two points, as illustrated in Figure 2. 
Thus, to identify candidate arcs, we divide the dual-space for the events (which 
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is also a sphere) into equal sized and shaped bins. A coarse example of the 
binning scheme [4] is shown in Figure 3. For our data analysis, we loop over 
all events in our data set. For each event, we increment those bins through 
which the corresponding great circle passes. The bin values therefore represent 
the total number of great circles that lie within its boundaries. The arcs are 
then identified by an excess in a particular bin. 

To match the dual-space statistics corresponding to the 271 events in our data 
set [3], we chose a 6° x 6° bin-size, giving a total of 4584 bins. For display pur- 
poses (figures 5, 6) the bins are numbered in latitude from north pole to south 
and around the axis. The dual-space bin contents are then compared to expec- 
tations from an entirely random, isotropically generated arrival distribution 
(without point or extended sources) to identify significant excesses. 
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Fig. 1. Simulated "polarization" plot by Gustavo Medina- Tanco [1]. It shows cosmic 
rays with energies between 10^^'^ and 10^'^ eV from point source. The events are 
deflected while passing through the galactic magnetic field forming arc shapes. While 
the events originate at point sources, the magnetic field deflects the events to form 
the arc shapes which one sees in the plot. 
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Fig. 2. Illustration of the intersection of five dual-space circles associated with five 
events that lie on an arc. The arrows point to the two intersection points. 




Fig. 3. An illustration of the sky divided into coarse, equal-sized and equal-shaped 
bins. A finer bin size is used in the actual analysis [4]. 



2 Applying The Hough Transform. 



We apply the Hough transform method to the stereo HiRes data to search 
for candidate arcs of events. The data selected are those events with energies 



4 



90" 





Fig. 4. The 271 event stereo data set in Hammer- Aitoff projection. 

> lO^^eV. Particles with energies < 10^^ eV are expected to be deflected by 
larger angles such that any structure would likely be lost. We exclude them 
as not to dilute any potential signals in the data set. This data set contains 
271 events collected between December 1999 and 31 January 2004. The event 
selection and quality cuts are identical to previous published study using the 
HiRes stereo data. Details can be found in [3] and [5]. 

Figure 4 shows the particle's arrival directions (in equatorial coordinates) for 
the 271 events used in this analysis in Hammer-Aitoff projection. After apply- 
ing the Hough Transform to the data set, the bin occupancy in dual-space is 
calculated. The resulting histogram is shown in Figure 5. 



3 Search for significant excess. 



In order to identify arcs, we look for apparent excesses over what is expected 
from an isotropic background. To estimate the significance of such excesses in 
each bin, we compared the data count to the distribution obtained from 2 x 10'' 
isotropic simulated data sets. Each simulated data set contains 271 events, and 
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Fig. 5. The bin occupancy values of normal great circles to the stereo data set . The 

sky was divided into 4584 bins of equal size (6° x 6° ) with the north pole at the 
zenith. The bins are arbitrarily numbered in latitude from north pole to south and 
around the axis. 



is based on a full detector simulation to account for the geometrical acceptance. 
In addition, the timing information is convolved to reproduce accurately the 
detector exposure in different directions of the sky over the detector hve-time. 
The simulation used in this study is described in more details in [3]. 

From the Monte Carlo sets, we calculate the probability for each bin to have 
a simulated occupancy greater than or equal to that found in the data. We 
define the significance, S, for a bin to be given by S = log(-p) and P = 
where N is the total number of simulated data sets, and M is the number of 
MC sets which recorded an occupancy greater or equal to that seen in the data 
bin. The resulting significance values are shown in Figure 6, where each bin 
is indexed by both a longitudinal (RA) and a latitudinal (DEC) bin number. 
Bins with S > 2.5 are designated as "high significance". 

For the seven high-significance bins we show the occupancy histogram for the 
2 X 10^ isotropic MC set in Figure 7, plotted in log scale. A vertical fine is drawn 
on each plot at the bin occupancy level found in the data. The probability of 
randomly finding at least this occupancy level in an isotropic data set is then 
the ratio of the area under curve to the right of the line to the whole area 
(2 X 10^). The significance "S" of each of the seven bins is shown in the upper 
right hand corner of each histogram. 
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Fig. 6. The significance values for each bin. As in figure 5, the bins are indexed by 
both a longitudinal (RA) and a latitudinal (DEC) bin number. 

4 Global Significance of data excess 



The next study we made was to determine whether the arcs found in the 
data are consistent with random fluctuations over the full sky, not just in the 
location in which they occurred. To do this, we compare our 2 x 10^ isotropic 
simulated sets to an independent set of 1000 isotropically simulated data sets, 
each again with 271 events. We applied the same analysis and cuts as we did 
to the real stereo data set. Including the "significance cut" of S> 2.5. We 
found that 99.6% of the time, the maximum significance of the isotropically 
generated data sets exceeded 2.5. Also, in 41.2% of the samples, the simulated 
data sets recorded a significance in at least one bin of > 3.5, which is the 
highest significance found in the real data set. Finally, the greatest number of 
events found on an arc in the real data is 19, as shown in Figure 7. This was 
also found to be quite common in the simulation, since arcs with S> 2.5 and 
at least 19 events were found in 90.4% of the simulated isotropic sets. 



5 Conclusion 



The Hough transform was applied to the HiRes stereo data and seven arcs were 
found with significance S > 2.5. However, when a global significance analysis 
was performed, similar results were commonly reproduced within isotropic 
data sets of equivalent size. Hence, the observation of seven arcs in the data 
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Fig. 7. Histogram of bin x from 2 x 10'^ simulated data sets for the seven most 
significant (S > 2.5) bins in the data. The corresponding data bin value is indicated 
in each case by a vertical line. The S values are shown in the upper right hand 
corner. 

is not significant. We conclude that we do not see any evidence of arc-like 
structures in the HiRes stereo data set. 
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